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SECTION  I 


INTRODUCTION 


Engine  components  have  traditionally  been  designed  using  a  crack  initiation  criterion. 
This  approach  has  been  very  successful  from  a  safety  standpoint,  but  the  conservatism  inherent 
in  this  method  has  resulted  in  poor  utilization  of  the  intrinsic  life  of  the  component.  The 
development  of  high  temperature  fracture  mechanics  has  permitted  basing  residual  life  analyses 
on  a  crack  propagation  criterion.  Using  this  approach,  lives  from  a  specified,  defect  size  can  be 
calculated  and  combined  with  periodic  inspection  to  determine  component  retirement. 

Propagation  analyses  classically  employed  in  residual  life  predictions  are  deterministically 
based.  They  typically  account  for  materials  scatter  by  the  incorporation  of  a  safety  factor.  A 
more  rigorous  treatment  of  materials  scatter  is  necessary  to  permit  maximum  utilization  of 
component  life. 

The  objective  of  this  program  is  to  develop  a  methodology  capable  of  incorporating  the 
scatter  observed  in  crack  growth  data  into  component  residual  life  analyses.  Such  a  methodology 
is  desired  for  use  in  damage-tolerant  and/or  Retirement  for  Cause  (RFC)  concepts.  This 
objective  includes:  identifying  the  distribution  functions  which  best  describe  crack  growth  rate 
(da/dn)  behavior;  characterizing  fatigue-crack  propagation  (FCP)  controlling  parameters  as  to 
their  effects  on  crack  growth  rate  variability;  studying  correlations  between  this  variability  and 
propagation  life  distributions;  and  developing  a  generic  methodology  applicable  to  all  engine 
materials,  although  particular  program  emphasis  is  placed  on  IN100,  Waspaloy,  and  Ti  6-2-4-6. 


When  the  critical  components  of  gas  turbine  engines,  such  as  disks,  must  be  used  beyond 
the  fatigue  crack  initiation  stage,  an  accurate  prediction  of  crack  growth  damage  accumulation 
for  these  components  under  service  loading  spectra  becomes  very  important.  Test  results  under 
laboratory  conditions  of  the  crack  growth  rates  for  engine  materials  such  as  IN100  and 
Waspaloy,  exhibit  considerable  statistical  variability  [e.g..  Refs.  1-7].  Typical  test  results  of  the 
crack  growth  rate  da/dN  versus  the  stress  intensity  range  AK  are  shown  in  Figures  1  and  2  for 
1N100  under  two  different  test  conditions.  As  a  result,  for  a  rational  prediction  of  the  fatigue 
crack  growth  damage  accumulation  of  enjgine  components,  it  is  essential  to  take  into  account  its 
inherent  statistical  variability.  Unfortunately,  a  statistical  model  and  analysis  procedure  for  the 
crack  growth  damage  behavior  of  engine  materials  have  not  been  established  to  date. 

Numerous  studies  have  been  carried  out  in  the  past  decade  on  the  crack  growth  behavior  of 
aluminum  alloys,  steels,  and  titaniums  (e.g.,  Ref.  8).  These  studies  are  essentially  for  airframe 
materials. 

i  • 

In  general,  literature  dealing  with  the  statistical  analysis  of  fatigue  crack  growth  damage 
accumulation  can  be  classified  into  three  categories.  Papers  in  the  first  category  consider  the 
random  crack  propagation  resulting  from  random  fatigue  loading  spectra  or  random  service 
loads,  such  as  narrow-banded  or  wide-banded  Gaussian  random  processes.  Typical  literature  is 
given  in  Ref.  9-18.  Within  this  category,  another  effort  has  also  been  made  to  convert  random 
service  loads  into  simplified  deterministic  loading  spectra  that  will  result  in  an  equivalent 
fatigue  crack  growth  behavior,  e.g.,  Ref.  19-23,  thus  saving  a  significant  amount  of  computation¬ 
al  effort.  A  conversion  of  cycle-by-cycle  integration  for  the  crack  growth  into  a  flight-by-flight 
integration  is  a  typical  example  (Ref.  19-23). 

Papers  in  the  second  category  deal  with  statistical  crack  propagation  in  the  small  crack  size 
region  using  fractographical  results  in  order  to  establish  the  so-called  equivalent  initial  flaw  size. 
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o.g.,  References  24-33.  Literature  in  the  third  category  deals  with  statistical  crack  growth 
behavior  under  deterministic  laboratory  loading  spectra,  e.g.,  Kef.  33-45,  thus  reflecting  the 
inherent  crack  growth  variability  due  to  the  material  itself. 

Unlike  airframe  structures,  the  loading  environments  of  engine  components  are  rather 
predictable.  Hence  literature  in  the  third  category  is  relevant  to  the  present  project  and  has 
been  examined  carefully.  While  the  statistical  characterization  of  the  crack  growth  rate 
presented  in  Ref.  34  and  lift  is  very  reasonable,  the  methodology  is  not  applicable  to  the  present 
problem  because  our  sample  size  is  too  small  under  each  test  condition.  Likewise,  the  statistical 
distributions  of  the  crack  size  in  service  have  not  been  solved  satisfactorily  in  Ref.  34  and  35. 
The  statistical  models  and  analysis  methodologies  presented  in  Ref.  33  were  developed 
essentially  for  aluminum  alloys  in  which  a  simple  crack  growth  rate  equation  can  be  used.  The 
statistical  model  presented  in  Ref.  38-41  appears  to  be  superior  to  that  of  Ref.  36-37.  However,  a 
considerable  amount  of  effort  in  further  research  and  development  is  needed  before  it  can  be 
applied  to  engine  components  due  to  the  complex  loading  environments. 

The  following  problems  are  unique  to  engine  materials  and  environments:  (1)  the 
statistical  distribution  of  the  crack  growth  damage  accumulation  is  influenced  by  many 
parameters  such  as  temperature  T,  loading  frequency  i\  stress  ratio  R,  holding  time  Th,  etc., 
(2)  the  number  of  specimens  tested  under  a  single  environmental  condition  mentioned  in  (1)  is 
very  small,  and  (3)  the  homogeneous  data  base  does  not  exist,  since  each  specimen  was  usually 
tested  using  different  specimen  geometry,  maximum  toad,  initial  crack  size,  and  final  crack  size. 
Thus,  it  is  necessary  to  develop  new  and  simple  statistical  models  capable  of  dealing  with  the 
unique  features  associated  with  engine  environments  on  the  basis  of  the  fracture  mechanics 
approach. 

Two  fracture  mechanics-based  statistical  models  for  the  fatigue  crack  growth  damage 
accumulation  in  engine  materials  are  proposed  and  investigated.  These  statistical  models  are 
based  on  the  synergistic  fracture  mechanics  models,  i.e.,  hyperbolic  sine  crack  growth  rate 
function,  developed  hv  Pratt  &  Whitney  [Ref  1  ].  Test  results  of  IN-100,  a  superalloy  used  in  the 
F100  engine,  at  various  elevated  temperatures,  loading  frequencies,  stress  ratios,  etc.,  have  been 
compiled  and  analyzed  statistically,  The  statistical  distributions  of(l)  the  crack  growth  rate;  (2) 
the  propagation  life  to  reach  any  given  ctack  size;  and  (3)  the  crack  size  at  any  service  life,  have 
been  derived.  It  is  demonstrated  that  the  correlation  between  the  test  results  and  two  statistical 
models  ! ,  very  good. 
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SECTION  II 


LOGNORMAL  crack  growth  rate  model 


The  crack  growth  rate  equation  proposed  by  Pratt  &  Whitney  Aircraft  (Ref.  1)  is  given  by 

Y  -  C,sinh  |C,(X  +  C,)|  +  C4 

in  which 

Y  -  log ,  K  -  lo*4K  (2) 

where  AK  -  stress  intensity  range,  Ct  is  a  material  constant,  and  C5,  C3and  C4are  functions  of 
temperature  T,  loading  frequency  v,  stress  ratio  R  and  others. 

The  first  statistical  model  assumes  that  the  crack  growth  rate  Y  is  a  homogeneous 
Gaussian  random  process 

Y  -  C|8inh  |C,(X  +  C3) ]  +  C4  +  Z(X) 

(3) 

in  which  Z(X)  is  a  homogenous  Gaussian  random  process  with  zero  mean  and  a  power  spectral 
density  S^w). 

The  statistical  model  described  by  Eq.  3  is  very  general  and  versatile.  Two  extreme  cases  of 
the  Gaussian  random  process  Z(X)  are  described  briefly  in  the  following: 

On  one  extreme,  Z(X)  is  a  Gaussian  white  noise.  For  a  Gaussian  white  noise,  the  correlation 
between  Z(X,)  and  Z(Xa)  is  zero  as  long  as  X,  *XS,  and  the  autocorrelation  function  is  a  Dirac 
delta  function  of  X.  Physically,  the  Gaussian  white  noise  resembles  the  random  walk  model. 
Sample  functions  of  Z(X)  and  Y(X)  are  schematically  shown  in  Figure  3. 
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Figure  $.  Sketch  of  Sample  Functions  of  Stationary  Gaussian  White  Noise  Z(X)  and 
Lor  Crack  Growth  Rate  Y(X) 

Based  on  the  white  noise  model,  the  statistical  dispersion  of  the  number  of  load  cycles, 
N(a),  to  reach  a  given  crack  size  “a"  or  the  crack  size,  a(n),  at  any  number  of  cycles,  n  is  the 
smallest  in  the  category  of  the  general  Gaussian  random  process.  Hence,  the  Gaussian  white 
noise  model  results  in  an  unconaervative  prediction  for  either  N(a)  or  a(n).  Consequently,  such 
an  extreme  case  is  not  appropriate  for  the  engineering  analysis  and  design  because  of  its 
unconservative  nature.  Likewise  the  mathematical  solution  for  such  a  problem  is  quite  involved. 

On  the  other  extreme,  Z(X)  iB  completely  correlated  for  all  X  values,  i.e.,  Z(X)=Z  is  a 
Gaussian  random  variable.  Physically,  the  random  variable  model  indicates  that  a  specimen 
starting  with  a  faster  (or  slower)  crack  growth  rate  will  always  preserve  the  faster  (or  slower) 
crack  growth  rate  throughout  its  crack  propagation  life.  Sample  functions  are  schematically 
shown  in  Figure  4  to  illustrate  the  behavior  of  Z(X)  and  Y(X). 
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Figure  4.  Sketch  of  Sample  Functions  of  Gaussian  Random  Variable  Z(X)  and  Log 
Crack  Growth  Rate  Y(X) 

Baaed  on  the  random  variable  model,  the  statistical  dispersion  of  N(a)  or  a(n)  is  the  largest 
in  the  general  cIbbs  of  Gaussian  random  processes.  Thus,  the  statistical  prediction  for  the  crack 
growth  damage  accumulation  using  such  a  model  is  conservative.  From  the  mathematical 
standpoint  the  random  variable  model  is  the  simplest  model.  As  a  result,  the  random  variable 
model  appeare  to  be  useful  and  appropriate  for  the  engineering  analysis  and  design  purpose,  and 
it  will  be  presented  in  this  report. 

In  reality,  sample  functions  of  the  crack  growth  rate  may  lie  between  these  two  extreme 
cases  as  schematically  shown  in  Figure  5.  Investigation  of  the  general  case  is  currently  underway 
and  the  results  will  be  presented  in  the  next  report. 
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Figure  5.  Sketch  of  Sample  Functions  of  Stationary  Gaussian  Random  Process  Z(X) 
and  Log  Crack  Growth  Rate  Y(X) 


Based  on  the  random  variable  model,  Bq.  3  can  be  written  bb 

Y  -  C|8inh[C2(X+C:,))  +  C„  +  Z  (4) 

where  Z  ia  a  normal  random  variable  with  zero  mean  and  a  standard  deviation  at. 

It  follows  from  Eqs.  4  and  2  that  the  log  crack  growth  rate  Y  is  a  normal  random  variable 
and  the  crack  growth  rate  da/dn  is  a  lognormal  random  variable  defined  in  the  positive  domain. 
The  mean  value,  My.  and  the  standard  deviation,  oy,  of  Y  can  be  obtained  from  Eq.  4  as 

My  -  C|sinh(C‘j(X+Cn)l  +  0,  (5) 

(Tv  =  av  (6) 

Hence,  the  probability  density  "unction  of  the  log  crack  growth  rate  Y  is  given  by 
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in  which  Myund  <»yare  k'vo.i  !:y  Eqs.  6  and  6. 

It  is  observed  from  Eqs.  5-7  that  the  probability  density  function  of  Y  involves  four 
parameters,  i.e.,  C2,  C;),  C4and  <ry.  These  parameters  will  be  estimated  from  test  results  using  the 
method  of  maximum  likelihood. 

Let  (x,,y,)  for  i  ■  1,2, . n  be  a  set  of  observations  of  X  and  Y  under  a  particular  test 

condition  (y,R,T)  where  v  *  loading  frequency,  R  -  stress  ratio  and  T  *  temperature.  The 
likelihood  function  L(y,,yj . ,yn;  Cj,C3,C4,<ry)  is  given  by 


■  1  J  1  /  Vi— My,  V21 

L(yi, . yn;C2,C3,C4,<ry)  —  a  2  y  ey  )  j 


(8) 


in  which  it  follows  from  Eq.  5  that 

My,  -  Cisinh[C2(xi+C3)]  +  C4 

Taking  the  natural  logarithm  of  both  sides  of  Eq.  8,  one  obtains 
8nL(yt . yn;C2,C3,C4,<ry)  -  -n8n(\^»  <jy)  -  (2ffy2)"1  S  (yt  -  My,)2 


(9) 


(10) 


The  maximum  likelihood  estimates  of  C8,  C3,  and  C4  and  «ryare  determined  by  maximizing 
Eq.  10  as  follows: 


3C2 

-  0 

an 

68nL(yiM...tyn;C2,C3,C4l<7y) 

(12) 

dC\\ 

-  0 

d8nL(yi,....,yn;C2,C;^C4,ay) 

**  0 

(13) 

dC4 

08nL(y  ymC^C;^  ,o  y) 

-  0 

(14) 
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Substituting  Eq.  10  into  Eqs.  11-14,  one  obtains  the  following  non!'  lear  coupled  equations 


n 


22  (  (Xi+C;,)(yj-/iyi)cosh[C2(x|+C3)l  =  0 

(16) 

£  j  (yi-Myi)co8htC2(xi+C3)l  -  0 

(16) 

^  (yt-My,)  -  0 

(17) 

>  T  . 

(18) 

in  which  /iy|  a  function  of  C,2,  Ca,  C4  and  x,  (Eq.  9). 

Equations  16-17  do  not  involve  <ry  and  hence  they  can  be  used  to  solve  for  C2,  Ca  and  C4 
simultaneously.  After  C.2  C.(  and  C4  are  determined  from  Eqs.  16-17,  the  standard  deviation  ay 
can  easily  be  computed  from  Eq.  18. 

It  should  be  mentioned  that  Eqs.  16-17  are  identical  to  those  equations  derived  by  using 
the  method  of  least  squares  as  described  in  Ref.  1.  Therefore,  the  least  square  estimates  of  C2, 
C.,  and  C4  are  identical  to  those  obtained  using  the  method  of  maximum  likelihood  described 
above.  Thus  the  computer  program  developed  in  Ref.  1  for  the  numerical  solution  of  Eqs.  16-17 
can  be  used. 

Since  Y  is  a  normal  random  variable,  the  crack  growth  rate  G  =  da/dn  follows  the 
lognormal  distribution.  The  coefficient  of  variation  of  the  crack  growth  rate  G,  denoted  by  V,  is 
related  to  the  standard  deviation  <ryof  Y  through 

V  -  [0<"y«"H»a  -  l  ]'»  (19) 

The  distribution  function  Ft.(»j)  and  the  probability  density  function  fC!(v)  of  the  crack  growth 
rate  G  ~  da/U(.  .’re  given  by 


in  which  *  ((•)  is  the  standardized  normal  distribution  function 

*  W  =  J^e  **/2  dt  (22) 
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nnd  nv  i«  given  by  Eq.  5 

My  -  C,8inhlC2(X+C:,)l  +  C<  (23) 

whereas  Ca,  C3,  C4and  iry  are  estimated  from  Eqs.  15-18. 

Teat  results  for  the  log  crack  growth  rate  Y  =  log  (da/dn)  versus  the  log  stress  intensity 
range  X  «  log  AK  are  presented  in  Figure  6  as  discrete  points  for  five  specimens  in  the  teat 
condition  No.  1  (T  ■  1000°F,  v  =  10  cpm,  R  "  0.1).  All  the  crack  growth  rate  data  shown  in 
Figure  6  have  been  used  to  estimate  C2,  C,,  C4  and  ey  using  the  method  of  maximum  likelihood. 
The  results  are  presented  in  the  first  row  of  Table  1  and  plotted  in  Figure  6  as  a  solid  curve.  The 
maximum  likelihood  estimates  of  C3,  C3,  C4  and  <ry  are  shown  in  Table  1  for  different  teat 
conditions  (v,  R,  T).  Also  presented  in  Table  1  are  the  corresponding  coefficients  of  variation  V, 
Eq.  19,  of  the  crack  growth  rate  G  ■  da/dn. 
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TABLE  1.  MAXIMUM  LIKELIHOOD  ESTIMATE  OF  C2,  Cv  C4,  STANDARD  DEVI¬ 
ATION  AND  THE  COEFFICIENT  OF  VARIATION  V  OF  Da/Dn 


Test  Condition * 

C, 

c. 

C, 

ft 

V 

Corf,  of 
Variation 

No.  of  Data 
f'ninls 

1 

3.8033 

-1.5239 

-4.3563 

0.1673 

40.00 ‘V 

258 

2 

4.9323 

-1.4073 

-3.9895 

0.1692 

40.49  r; 

150 

3 

4.3093 

-1.3032 

-4.4450 

0.1240 

29.1417 

130 

4 

3.8624 

-1.5054 

-4.1594 

0.2078 

50.72 '7- 

188 

5 

3.8982 

-1.5376 

-3.9341 

0.1026 

23.96 '7 

342 

Average  36.86 'V 


*Test  Conditions 


1. 

T  =  1000-F, 

v  — 

10  cpm.  R  -  0.1 

2. 

T  -  1350°F, 

V  — 

10  cpm.  R  =  0.1 

3. 

T  =  lZOO^F, 

V  — 

10  cpm.  R  =  0.5 

4. 

T  =  1200°F. 

v  — 

20  cpm,  R  -  0.05 

5. 

T  =  1200°  F, 

V  — 

10  cpm.  R  =  0.1 

To  show  the  validity  of  the  assumption  that  Z  follows  the  normal  distribution  with  zero 
mean,  sample  values  of  Z  corresponding  to  test  results  (x^Vj),  denoted  bv  Z;,  are  computed  from 
Eq.  4  as 


z,  =  y*  -  Cjsinh  [G>(xi+C;,)]  -  C4;  for  i  =  1,2 . n 


(24) 


Sample  values  of  Zj  (i  =  1,2 . n)  for  each  test  condition  (v,  R,  T)  are  plotted  on  the  normal 

probability  paper  in  Figures  7-11  as  circles.  The  normal  probability  paper  is  constructed  using 
linear  scales  for  both  the  ordinate  Q  and  the  abscissa  Z.  Sample  values  of  Z  are  ranked  in  an 
ascending  order  and  the  plotting  position  for  the  jth  data  point  z]  is  Q3  =  4>"!  (j/n+1).  Also  shown 
in  these  figures  are  the  straight  lines  representing  the  normal  distribution  of  Z  with  zero  mean 
and  a7  =  iry  given  in  Table  1.  It  is  observed  in  each  figure  that  the  test  results  z.  are  scattered 
around  the  respective  straight  line,  indicating  that  the  normal  distribution  is  reasonable  for  Z. 
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Figure  II.  Normal  Probability  Plot  for  Z  ( Test  Condition  5) 


Tests  for  goodness-of-fit  using  the  Kolmogorov-Smirnov  (K-S)  test  have  been  carried  out. 
The  observed  K-S  statistics  Dn  along  with  the  critical  value  D„“  associated  with  a  level  of 
significance  are  presented  in  Table  2  for  each  test  condition  (t,  R,  T).  From  Table  2  the  normal 
distribution  is  acceptable  for  Z  at  a  20 'V  level  of  significance  for  test  conditions  3-5,  whereas  the 
level  of  significance  is  5'V  for  test  conditions  1  and  2. 


TABLE  2.  COMPARISON  OF  OBSERVED  Dn  AND  CRITICAL  VALUES  OF 
Dn"  IN  THE  KOLMOGOROV-SMIRNOV  TEST  FOR  Z 


Tent 

Condition 

Sample 

Site,  n 

Observed 

[)„ 

Crtilical  Values  of  D„" 

«  =  0.2 

«  -  0.1 

u  -  0.05 

a  -  0.01 

1 

258 

0.081 

0.067 

0.076 

0.086 

0.102 

2 

150 

0.114 

0.088 

0.100 

0.111 

0.133 

:t 

i:to 

0.091 

0.094 

0.107 

0.119 

0.143 

4 

188 

0.064 

0.078 

0.089 

0.099 

0.119 

fi 

342 

0.050 

0.068 

0.066 

0.736 

0.088 

n  -  level  of  HiicnificBni’e 
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It  is  important  to  emphasize  that  the  present  statistical  model  (random  variable  model  in 
Kq.  A),  referred  to  as  the  first  statistical  model,  has  a  significant  advantage  in  that  the  statistical 
distribution  of  the  crack  growth  damage  accumulation  can  be  obtained  analytically. 

Statistics!  Distribution  of  Crack  Growth  Rata 

After  C.j,  C.,,  C4and  <rv  have  been  estimated  from  Eqs.  15-18  for  each  test  condition  given  in 
Table  1,  the  distributions  of  the  crack  growth  rates  of  Y(AK,Ct)  =  log  da/dn  and  G(AK,Cj)  = 
da/dn,  for  i  =  2,3,4  are,  respectively,  normal  and  lognormal  as  given  by  Eqs.  7, 20-21. 

To  theoretically  demonstrate  the  behavior  of  the  distributions  of  Y  and  da/dn,  let  zy  be  the 
7  percentile  of  Z,  i.e., 

7%  =  P[Z>z,|  =1-4.  (25) 

where  Z  follows  the  normal  distribution  with  zero  mean  and  the  standard  deviation  nt  =  oy  (see 
Eq.  6).  Inversely  Eq.  25  can  be  written  as 

z,  =-  e*  (1-7%)  (26) 

in  which  <rt=<;yand  4>  1  (  )  is  the  inversed  function  of  4>  (  ). 

The  7  percentile  of  the  log  crack  growth  rate  yy  (AK.C,),  (i  =  2,3,4),  follows  from  Eq.  4  as 

y,  (AK.Cj)  =  (>inh|C,(X  +  C:,)|  +  C4  +  z,  (27) 

and  the  7  percentile  of  the  crack  growth  rate 


(28) 


becomes 

G,  (AK.Cj)  =  (1()).v,<^K.(M  (29) 

in  which  yy  (AK,C()  is  given  by  Eqs.  26  and  27. 

Note  that  G  (AK,C,)  and  y^AK.C,)  are  functions  of  the  stress  intensity  range  AK  and  C,  (i 
2, 3, 4)  that  depend  on  test  conditions  (i\  R,  T). 

The  7  percentile  log  crack  growth  rate  path  yy  (AK,  (’>)  are  computed  from  Eqs.  26  and  27 
and  plotted  in  Figure  12  for  the  test  condition  No,  1  (i>  10  cpm,  R  -  0.1  and  T  --  1000°F)  for  7 

5,  25,  50,  75,  95.  Curves  in  Figure  12  also  represent  possible  sample  paths  of  the  crack  growth 
rate.  For  instance,  the  crack  growth  rate  path  associated  with  7  -  5',  indicates  that  5'V  of  the 
specimens  will  have  a  crack  growth  rate  path  higher  (or  faster)  than  that  shown  by  the  curve.  It 
is  observed  from  Figure  12  that  sample  paths  of  the  crack  growth  rate  do  not  intermingle,  and 
that  a  specimen  with  a  higher  crack  growth  rate  always  preserves  it  over  the  entire  range  of  AK. 
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Figure  12.  Fcrcvntilcs  of  log  C 'rack  Growth  Hate  Y  as  Function  of  Lon  Stress 
Intensity  Ran, i>e  X  for  Model  I  and  Test  Condition  No.  I 


Statistical  attribution  of  Crack  Slxo  and  Tima  to  Roach  Qlvan  Crack  8ixa 

The  7  percentile  crack  size  after  n  load  cycles,  denoted  by  a^n),  can  be  computed  using  the 
7  percentile  crack  growth  rate  G^AK.C,)  as  follows: 

( ■&),-  wo) 

Aaj(7)  -  G^AK.CMnj  (31) 


m 

a>(n)  =■  ao+  2:  Aaj(7) 


(32) 


in  which  a()is  the  initial  crack  size  and 


The  cycle-by-cycle  numerical  integration  given  by  Eqs.  30-33  is  deterministic  and 
straightforward.  Hence,  by  varying  the  value  of  the  y  percentile,  one  obtains  from  a  cycle-by¬ 
cycle  integration  a  set  of  crack  growth  curves  a^n),  i.e.,  the  crack  size  versus  the  number  of 
cycles  for  each  y  value.  The  y  percentiles  of  the  crack  growth  damage  accumulation  a^n)  are 
plotted  in  Figure  13  for  the  test  condition  No.  1  (i*  =  10  cpm,  R  -  0.1  and  T  =  1000°F)  with  the 
following  specimen  geometry  and  loading  condition:  ASTM  CT  specimen,  thickness  =  0.25  in., 
width  =  2.602  in,,  maximum  load  =  1.4  ksi,  initial  crack  size  =  0.5  in.,  final  crack  size  2.0  in.,  C,  = 
0.5  and  C,,  i=2,3,4  given  in  Table  1.  For  instance,  the  crack  growth  curve  a7(n)  corresponding  to 
y -10  in  Figure  13  indicates  that  there  are  10',  of  the  specimens  having  a  crack  size  larger  than 
those  indicated  by  this  curve  at  any  number  of  load  cycles  n. 

TEST  CONDITION  NO.  1 


NUMBER  OF  CVCLES,  10**4  «*»«> 

Figure  Vi.  Distribution  for  Crack  Growth  Damage  Accumulation  a(n)  for  Model  I  and 
Test  Condition  No.  / 


After  constructing  a  series  of  crack  growth  damage  accumulation  curves  a  (n)  for  many 
values  of  7  as  shown  in  Figure  13,  one  can  establish  (i)  the  distribution  function  FN(1,  (n)  of  the 
number  of  load  cycles  N(a)  to  reach  any  crack  size  “a”  by  drawing  a  horizontal  line  through  "a” 
and  (ii)  the  distribution  function  F,tnl(x)  of  the  crack  size  a(n)  at  any  number  of  load  cycles  n  by 
drawing  a  vertical  line  through  a. 

For  instance,  the  distribution  FN(<|  (n)  for  the  number  of  load  cycles  to  reach  the  crock  sizes 
u  1.0  in.  and  2.0  in.  are  shown  in  Figure  14(a)  and  14(b),  respectively,  as  solid  curves.  The 
complement  of  the  distribution  function  of  the  crack  size  a(n)  at  n  25,000  load  cycles,  i.e., 
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F*B|11>  (x)  1  K<|fl|  (x),  is  presented  in  Figure  15  as  a  solid  curve.  In  Figure  15  the  ordinate  F*,(n)is 
the  probability  that  the  crack  size  will  exceed  the  corresponding  value  indicated  by  the  abscissa. 
For  instance,  the  probability  is  0.41  that  the  crack  size  at  n  -  26,000  cycles  will  exceed  0.65  in.  as 
shown  by  a  square  on  the  solid  curve  in  Figure  15.  Figure  15  is  referred  to  as  the  crack 
exceedance  curve.  In  conclusion,  on  the  basis  of  the  present  statistical  model,  the  distribution  of 
the  crack  growth  damage  accumulation  presented  in  Figure  15  contains  all  the  information 
needed  for  the  analysis  of  retirement  for  cause  of  engine  materials. 

Correlation  With  Tatt  Raaulta 

One  of  the  difficult  problems  encountered  in  the  present  program  is  that  the  number  of 
test  specimens  under  each  environmental  condition  (v,  R,  T)  is  very  small.  Even  with  such  a 
small  sample  size,  each  specimen  is  usually  tested  using  different  specimen  geometry,  initial 
crack  size,  final  crack  size  and  maximum  load.  A  summary  of  specimen  geometries  and  loading 
conditions  for  all  the  test  specimens  is  presented  in  Table  3.  It  is  obviouB  from  Table  3  that  the 
experimental  results  of  the  crack  growth  damage  accumulation  are  not  homogeneous,  making 
the  correlation  study  difficult.  A  qualitative  correlation  study  is  carried  out  in  the  following. 

In  order  to  correlate  the  test  results  with  the  predictions  based  on  the  statistical  model 
proposed  previously,  homogeneous  test  environments  have  been  assumed  for  each  test  condition 
as  shown  in  Table  4.  The  maximum  loads  PmH  given  in  Table  4  are  chosen  in  order  to  avoid 
excessive  extrapolation  far  into  the  region  of  AK  in  which  actual  test  results  do  not  exist.  Test 
results  of  the  crack  growth  rate  for  each  specimen  are  best-fitted  by  Eq.  1  using  the  method  of 
least  squares  to  estimate  values  of  Cv  C.,and  C4.  The  least  squares  beBt  fit  procedures  have  been 
described  in  Ref.  1.  Sample  values  of  C,  (1 =2,3,4)  thus  obtained  for  each  specimen  are  presented 
in  Table  6  for  various  test  conditions.  As  an  example,  the  best-fitted  crack  growth  rate  for  the 
five  specimens  in  the  test  condition  No.  1  are  plotted  in  Figure  16  using  Eq.  1  and  C,  0=2,3, 4) 
values  given  in  Table  5.  The  regions  of  Y  and  X  shown  by  solid  curves  are  the  regions  in  which 
test  data  exist,  whereas  the  dotted  curves  represent  the  extrapolation  into  the  regions  where  test 
data  are  nonexistent, 
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DISTRIBUTION  FUNCTION  DISTRIBUTION  FUNCTION 


T.  C.  1 


fA 


/O  a  =  1.0  INCH 

- MODEL  2 

- MODEL  1 

O EXTRAPOLATED 

TEST  RESULTS 

— i 1 1 1 J I I l 

4  8  12  16  20 

NUMBER  OF  CYCLES,  104 


T.  C.  1 


a*  2.0  INCH 
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Figure  N  Distributions  of  Cycles  to  Reach  a  Given  Crack  Size  for  Test  Condition 


PROBABILITY  OF 
CRACK  EXCEEDANCE 


FD  235872 


Figure  15.  Probability  of  Crack  Exceedance  at  25,000  Load  Cycles  for  Model  1  and 
Test  Condition  No.  1 
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TABLE  3.  SPECIMEN  GEOMETRY  AND  MAXIMUM 
LOAD  FOR  EACH  TEST  SPECIMEN 


Specimen 

No. 

fl  «» 

(in.) 

ar 

(in.) 

IV 

(in.) 

H 

fin.l 

f>™, 
f kip a) 

Tent  Condition  No.  1 

0582 

0.812 

1.788 

2.502 

0.250 

1.008 

1082 

0.769 

1.819 

2.600 

0.499 

1 ,6411 

10811 

0.742 

2.021 

2.6118 

0.501 

1.8211 

1828 

0.724 

1.608 

2.1X14 

0,497 

1.846 

1820 

0.1170 

1.427 

2.008 

0.499 

3.664 

Test  Condition  No.  2 

01118 

0.8188 

1.11711 

2.497 

0.479 

1.706 

01171 

0.838 

1.797 

2.409 

0,502 

2.027 

0541) 

0.804 

1.757 

2.498 

0,462 

1,754 

1705 

0,617 

1,840 

2.6011 

0.601 

2.5711 

Teat  Condition  No.  3 

011114 

0.8110 

1.0110 

2.408 

0.501 

2.0 18 

011011 

0,028 

1.027 

2.606 

0.600 

2.037 

1826  0.020  1.407  2.00  0.500  2.20 


1827 

0.980 

1.618 

2,001 

0.500 

1.628 

Test  t'ondition  No,  4 

0710 

0.860 

1.21*1 

2.400 

0.828 

11.146 

0711 

0,076 

1.710 

2.513 

0.400 

1.076 

0717 

1,168 

1.911 

2.500 

0.364 

0.916 

11186 

0.026 

1.948 

2,507 

0.408 

2.210 

11189 

0.818 

1.805 

2.506 

0.600 

2.426 

1722 

0.162 

0.545 

1.261 

0.11)8 

2.4118 

Test  Condition  No.  S 

05114 

0.896 

2.012 

2.605 

0,501 

1.601 

0641 

0,808 

1.717 

2,406 

0.409 

2.6011 

OfllillAK 

1.006 

1.71111 

2.512 

0.260 

O.lXXl 

0686 AK 

0.828 

1.7(18 

2.507 

0.870 

11.620 

10511 

0.862 

1.76(1 

2.403 

0.290 

1.480 

11141 

0.088 

0.11511 

0.986 

0.29(1 

7.168 

0114(1 

1.401 

2.1X10 

2.6011 

0.501 

1.1156 

10711 

0.461 

1.801 

2,407 

0.601 

11,6311 

1674 

0.476 

1.880 

2.601) 

0.400 

3.811 

TABLE  4. 

ASSUMED  HOMOGENOUS  TEST  ENVIRON- 

MENTS  FOR  TEST 

SPECIMENS 

Tent 

a  «)) 

<V 

IV 

H 

/' 

turn 

Condition 

(in.) 

(in.) 

fin.> 

<  in.) 

(kips) 

1 

0.5 

2.0 

2.602 

0.25 

1.4 

2 

0.5 

2.(1 

2.5011 

0.601 

2.6 

11 

0.5 

1.8 

2,00 

0.50 

2.00 

4 

0.6 

2.(1 

2.500 

0.50 

11.0 

6 

0.6 

2.(1 

2  497 

(1.601 

11.51111 

aldl  Initial  crack  Hire 

ty  -  final  crack  hIm* 

Hpeciman  width 
upeclmen  thicknemt 
I* maaimum  load 


TABLE  r>. 


TEST  RESULTS  OP  Cr  C,  AND  C4  FOR 
EACH  SPECIMEN  UNDER  VARIOUS 
TEST  CONDITIONS 


Specimen 

No.  C.j 

c, 

(a) 

T*»t  Condition  No.  1 

T  - 

KMWK, 

V  “ 

10  cpm,  K  ”  0.1 

0582 

1.1820 

-1.145 

-4,690 

1082 

5.051 

-1.158 

-4.676 

108,1 

1.989 

-1.490 

-4,490 

1828 

1,962 

-1.517 

-4.146 

1829 

4.041 

-1.605 

-4.148 

tb> 

Tent  Condition  No.  2 

T  - 

1150* K. 

I.  - 

10  cpm,  R  »■  0.1 

01118 

5,149 

-1,097 

-4.810 

01171 

4.827 

-1.188 

-4.085 

05411 

1.821 

-1.518 

-1.741 

1705 

5.556 

-1.407 

-1.918 

(t) 

Tent  Condition  No.  0 

t  - 

12(H)'K, 

ii  “ 

10  cpm,  R  -  0,5 

0.1.14 

4.800 

-1.158 

-4.870 

oi«:i 

5.214 

-1.202 

-4.081 

182(1 

4.016 

-1,147 

-4.110 

1827 

4,704 

-1.280 

-4.511 

(d) 

Tent  Condition  No.  4 

T  - 

HOOT, 

,.  - 

20  cpm,  K  -  0,05 

0710 

5.554 

1.411 

4.276 

0711 

4.470 

-1.107 

4.116 

0717 

4.201 

1.501 

1.904 

1 185 

4.484 

1.560 

-4,140 

1189 

1.711 

4.400 

4.144 

1722 

4.110 

- !  .107 

•4.100 

(«> 

Tent  Condition  No.  5 

T  - 

1200*  K. 

v  « 

10  cpm,  K  ■»  0.1 

0514 

4.272 

-1.528 

-1.825 

0541 

1.780 

1.500 

1.858 

OOllAK 

1.122 

1.049 

1.711 

0085AK 

1.700 

1.500 

-  1,877 

1051 

1.804 

1.519 

-1.971 

H4I 

4.047 

1,511 

1,972 

0.14(1 

1.674 

1.109 

4.402 

1071 

1.101 

1.497 

1.960 

1074 

4.120 

1.510 

4,045 
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Figure  i 6',  Best-Fitted  Crack  Growth  Kate  for  each  Specimen  in  Test  Condition  No.  I 


The  best-fitted  C,,and  (Rvalues  lor  each  teat  specimen  presented  in  Table  5  are  used  to 
grow  the  crack  under  assumed  homogeneous  test  environments  shown  in  Table  4.  The  crack 
growth  damage  accumulation  thus  obtained  is  referred  to  as  the  “extrapolated  test  results"  in 
the  sense  that  they  are  not  the  actual  test  results. 

Let  us  consider  the  test  condition  No.  1.  Based  on  the  theoretical  model,  the  maximum 
likelihood  estimates  of  C2,  C.t,  C4  and  <ry  obtained  in  Table  1  are  used  to  construct  various 
percentiles  of  the  log  crack  growth  rate  Y  =  log  (da/dn)  versus  the  log  stress  intensity  range  X  = 
log  AK.  The  results  are  presented  in  Figure  17(a).  The  least  square  estimates  of  C2,  C.,and  C4for 
each  specimen  obtained  in  Table  5  are  used  along  with  Eq.  1  to  construct  Y  versus  X  as  shown  in 
Figure  17(b).  Figure  17(b)  is  referred  to  as  the  “extrapolated  test  results”  since  the  actual  test 
data  do  not  exist  in  the  entire  region  of  AK, 

With  the  assumed  homogeneous  test  conditions  given  in  Table  4,  the  log  crack  growth  rates 
displayed  in  Figure  17  have  been  integrated  to  yield  the  crack  growth  damage  accumulation. 
The  corresponding  crack  size,  a(n),  versus  the  number  of  load  cycles,  n,  is  shown  in  Figure  18.  A 
comparison  between  Figures  18(a)  and  18(b)  indicates  that  the  correlation  between  the 
extrapolated  test  results  and  the  statistical  model  is  reasonable. 

In  the  same  manner,  the  log  crack  growth  rate,  Y,  versus  the  log  stress  intensity  range,  X, 
and  the  crack  size,  a(n),  versus  the  number  of  load  cycles,  n,  have  been  computed  and  presented 
in  Figures  19-26  for  test  conditions  No.  2  to  5.  it  is  observed  from  these  figures  that  the 
correlation  between  the  extrapolated  test  results  and  the  statistical  model  is  reasonable. 

The  distributions  FN(|1)  (n)  for  the  number  of  load  cycles  N(a)  to  reach  any  crack  size  based 
on  the  statistical  model  is  plotted  in  Figures  27-31  as  solid  curves  for  all  the  test  conditions  and 
two  particular  crack  sizes.  The  corresponding  distribution  functions  for  the  extrapolated  test 
results  are  determined  from  Figures  18(b),  21{b)-24(b)  and  shown  in  Figures  27-31  as  circles.  In 
construct  ing  the  distribution  function  (circles)  of  the  extrapolated  test  results  in  Figure  27,  for 
instance,  one  obtains  five  data  points  from  Figure  18(b)  by  drawing  a  horizontal  line  through  a  ■= 
1.0  in.  These  five  data  points  are  ranked  in  an  ascending  order  and  the  ordinate  of  the  ith  data 
point,  is  i/(m+T)  where  m~-5  is  the  total  number  of  data  points.  Figures  27-31  demonstrate  a 
good  correlation  between  the  statistical  model  (solid  curves)  and  the  extrapolated  test  results 
( ■  ircles). 

The  procedures  for  determining  the  distribution  of  the  crack  size  a(n)  at  any  number  of 
load  cycles,  n,  have  been  described  previously.  The  complement  of  the  distribution  function 
F*„n(  (x)  of  n(n),  referred  to  as  the  crack  exceedance  curve,  is  plotted  in  Figures  32-36  as  solid 
curves.  For  instance,  the  square  shown  in  Figure  32  indicates  that  at  na 25,000  load  cycles  the 
probability  that  the  crack  size  will  exceed  0.65  in.  (abscissa)  is  0,41  (ordinate)  for  the  test 
condition  No.  1.  It  is  observed  that  the  exceedance  probability  reduces  as  the  crack  size 
increases.  The  extrapolated  test  results  for  a(n)  are  obtained  from  Figures  18(b),  23(b)-26(b)  by 
drawing  a  vertical  line  through  the  number  of  load  cycles,  n,  of  interest.  They  are  plotted  in 
Figures  32-36  as  circles  for  all  test  conditions.  Figures  32-36  demonstrate  that,  the  correlation  for 
the  crack  size  distribution  between  the  statistical  model  and  the  extrapolated  test  results  is  very 
reasonable. 
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h Hurt'  18.  (.'rack  Urowth  Damage  Accumulation  for  Test  Condition  No.  /,'  (a) 
Statistical  Model  and  (b)  Extrapolated  Test  Results 
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Fiuure  21.  Log  Crack  Growth  Rate  Y(X)  vs  Lag  Stress  Intensity  Range 
Condition  No.  4;  (a)  Statistical  Model  No.  1  and  (b)  Extrapt 
Results 
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b inure  22.  Lug  Crack  Growth  Rate  Y(X)  us  Log  Stress  Intensity  Range  X  for  Test 
Condition  No.  5;  (a)  Statistical  Model  No.  I  and  (b)  Extrapolated  Test 
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Figure  23.  Crack  Growth  Damage  Accumulation  for  Test  Condition  No.  2;  (a) 
Statistical  Model  No.  I  and  (b)  Extrapolated  Test  Results 
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Figure  24.  Crack  Growth  Damage  Accumulation  for  Test  Condition  No.  ,'J;  fa) 

Statistical  Model  No.  I  and  (b)  Extrapolated  Test  Results 
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Figure  25.  Crack  Growth  Damage  Accumulation  for  Test  Condition  No.  4;  (a) 

Statistical  Model  No.  I  and  (b)  Extrapolated  Test  Results 
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DISTRIBUTION  FUNCTION  DISTRIBUTION  FUNCTION 


NUMBER  OF  CYCLES,  104 


Figure  27.  Distribution  of  Cycles  to  Reach  a  (liven  Crack  Site  for  Test  Condition  No.  1 
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DISTRIBUTION  FUNCTION  DISTRIBUTION  CYCLES 
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Fifture  2H.  Distribution  of  Cycles  to  Reach  a  ( liven  Crack  Size  for  Test  Condition  No,  2 
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b ifiurv  29.  Distribution  of  Cycles  to  Reach  a  (liven  Crack  Size  for  Test  Condition  No,  3 
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DISTRIBUTION  FUNCTION  DISTRIBUTION  FUNCTION 


hi  Hurt'  HO.  Distribution  of  Cycles  to  Reach  a  Given  Crack  Size  for  Test  Condition  No. 
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Figure  ID,  Distribution  of  Cycles  to  Reach  a  (liven  Crack  Site  for  Test  Condition  No  5 
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Figure  ,‘12.  Probability  of  Crank  Exceedance  at  n~ 25,000  Load  Cycles  far  Test 
Condition  No,  1 
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Figure  bit.  Probability  of  Crack  Exceedance  at  n-  10,000  Load  Cycles  for  Test 
Condition  No.  2 
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PROBABILITY  OF  PROBABILITY  OF 

CRACK  EXCEEDANCE  3  CRACK  EXCEEDANCE 
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lure  34.  Probability  of  Crack  Exceedance  at  n^ 40,000  Load  Cycles  for  Test 
Condition  No.  3 
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Figure  35.  Probability  of  Crach  Exceedance  at  n~  10,000  Load  Cycles  for  Test 
Condition  No.  4 


SECTION  III 


RANDOMIZATION  OF  PARAMETERS 


Statistical  Modal 


The  second  statistical  model  is  expressed  as 


Y  -  C,sinh|C2(X  +  C;,)|  +  C4 


(34) 


in  which  Ca,  C.,  and  C4  are  considered  as  random  variables,  i.e.,  values  of  C2,  C3,  and  C4  vary  from 
specimen  to  specimen.  Thus  each  test  specimen  results  in  a  set  of  sample  values  of  C2,  C,  and  C4. 

To  characterize  the  statistical  distributions  of  C2,  C;1  and  C4  as  well  as  their  correlations, 
test  results  of  crack  growth  rates  for  each  specimen  are  best  fitted  by  Eq.  34  using  the  method  of 
least  squares  to  estimate  sample  values  of  C2,  C3  and  C4.  The  least  squares  best  fit  procedures 
have  been  discussed  in  Ref.  1.  Sample  values  of  C,  (i= 2,3,4)  are  approximated  by  their  estimates 
and  are  presented  in  Table  5  for  various  test  conditions.  For  instance,  the  best-fitted  crack 
growth  rate  for  five  specimens  in  the  test  condition  No.  1  are  plotted  in  Figure  16  UBing  Eq.  34  as 
well  as  the  C,  (i=2,3,4)  values  given  in  Table  5. 

The  mean  values  (p2,g.,,p4),  the  standard  deviations  (ff2,<r3,<74),  and  the  coefficients  of 
variations  (V2,V3,V4)  of  those  data  given  in  Table  6  have  been  computed  in  the  following 
manner: 
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V,  =  ffj/Mj,  V;,  =  V.|  =>  ff4/M4 


(37) 


in  which  (C21,C3l,C4i)  are  the  ith  sample  values  of  (C2C3,C4)  associated  with  the  ith  specimen,  and 
m  is  the  total  number  of  samples  (specimens)  for  a  test  condition. 
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The  correlation  coefficients  pm,  pM,  p42  among  Cj  (j- 2,3,4)  are  obtained  as  follows: 
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(39) 

(C4t~ M4>(C2i— Ma) 

(40) 

values  of  p2a,  pM  and  p42  thus  computed  are  also  given  in  Table  6  for  various  test  conditions.  To 
show  the  correlations  among  C2,  Cn  and  C4  under  different  test  conditions,  test  data  presented  in 
Table  5  are  shown  in  Figure  37  for  C2vs  C3,  C.,  vs  C4and  C4vs  C2. 


TABLE  6.  MEAN  VALUES,  STANDARD  DEVIA¬ 
TIONS,  COEFFICIENTS  OF  VARIATION 
AND  CORRELATION  COEFFICIENTS  OF 
C2,  C,  AND  C4 


Test 

Condition 

c. 

Mean 

Values 

1 

4.081 

-1,467 

-4.470 

2 

4.88ft 

-1.1168 

-4.161 

:t 

4.718 

-1.247 

-4.600 

4 

4.480 

-1.469 

-4,215 

6 

11.802 

-1.629 

-3.960 

.Standard 

Deviations 

1 

0.544 

0.101 

0.206 

2 

0.670 

0.181 

0.414 

11 

0.4112 

0.0711 

0.208 

4 

O.fidl 

0.092 

0.134 

ft 

0.299 

0.070 

0.180 

Coefficient* 

0/  Variation 

1 

0.11111 

-0.069 

-0.048 

2 

0.1117 

-0.119 

-0.100 

:i 

0.091 

-0.058 

-0.045 

4 

0.127 

-0.061) 

-0.01)2 

5 

0.079 

-0.048 

-0.046 

Correlation 

Coefficients 

Test 

Condition 

I’jj 

P. 11 

1‘ij 

1 

0.1211 

-0.982 

-0.121 

2 

0.651 

-  0.992 

-0.662 

:i 

0.7711 

-0.990 

-0.721 

4 

0.088 

0.918 

-0.264 

r» 

0.162 

-0.91)8 

-0.047 

From  the  results  preKented  in  Table  6  and  Figure  .'17  it  in  clear  that  the  number  of 
specimens  tested  under  each  condition  (t,  R,  T)  is  too  small  to  characterize  the  statistical 
distributions  of  C2,  C.,  and  C4.  Nevertheless,  from  the  analysis  results  presented  above,  several 
interesting  trends  have  been  observed  as  follows: 

(1)  The  mean  values  (p2,p3,p4)  are  functions  of  the  test  condition  (v,R,T). 

(2)  The  statistical  dispersions  or  the  coefficients  of  variation  (V2,V3,V4)  are  fairly  con¬ 
stant  over  various  test  conditions,  such  that  (V2,V3,V4)  may  be  assumed  to  be  independent  of  the 
test  condition  (v.R.T). 

(3)  The  correlations  p23  and  p4ii  between  C2  and  C3,  and  C2  and  C4  are  very  small  except  in 
the  third  test  condition,  indicating  that  C2is  almost  uncorrelated  with  C3and  C4. 

(4)  The  correlation  p34  between  C.,  and  C4  is  almost  unity,  indicating  that  C3  and  C4  may 
be  linearly  related. 

Because  of  the  difficulty  associated  with  the  small  sample  size  as  well  as  the  observations 
made  above,  the  following  assumptions  are  made  so  that  all  the  test  results  under  different  test 
conditions  can  be  pooled  together  to  characterize  the  distributions  of  C2,  C3and  C4. 

(1)  The  mean  values  (/i2,p3,p4)  of  C2,  C3and  C4  are  functions  of  (i\R,T),  i.e„  (p2,p3,p4)  vary 
from  test  condition  to  test  condition  and  their  values  are  given  in  Table  6. 

(2)  The  coefficients  of  variation  (V2,V3,V4)  of  C2,  C3  and  C4  are  independent  of  (u,R,T), 
i.e.,  (V2,V3,V4)  are  constants. 

Based  on  these  assumptions  we  normalize  all  the  Bample  values  of  C2,  C3  and  C4  by  their 
respective  mean  values  (p2,/i.,,p4).  The  resulting  data  will  have  the  mean  value  equal  to  unity  and 
the  same  standard  deviation.  Thus,  all  the  normalized  data  under  different  test  conditions  can 
be  pooled  together  to  determine  (1)  the  coefficients  of  variation  (V2,V,,V4),  (2)  the  correlation 
coefficients  among  C2,  C.,and  C4,  and  (3)  the  statistical  distributions  of  C2,  C,and  C4. 

Let  nij  be  the  total  number  of  test  specimens  under  the  jth  test  condition  (t»j,Rj,Tj)  and  J  be 
the  total  number  of  test  conditions.  Then  the  total  number  of  test  specimens  M  is  given  by 

.1 

M  =  S  m,  (41) 

.i  i 


Let  C2,  C.,and  C4  be,  respectively,  the  normalized  random  variables  of  C2(u,R,T),  C3(u,R,T) 
and  C4(u,R,T)  with  respect  to  their  mean  values,  i.e., 

C,  =  C-(u,R,T)/p,,(i»,R,T) 

0;t  -  C;|(u,R,T)/p;|(l\R,T)  (42) 

C.|  =  Cifu.R.Tl/pilv.R.T) 
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Test  Condition  No.  3 


Teat  Condition  No.  4 


Test  Condition  No.  5 
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Figure  37c.  Correlation  Between  C3  and  C4  for  Various  Test  Conditions 


Sample  values  of  C2,  C,  and  C4  are  obtained  by  dividing  each  data  point  in  Table  5  by  its 
corresponding  mean  value.  The  resulting  data  are  pooled  together  and  presented  in  Table  7.  The 
coefficients  of  variation  (V2,V3,V4)  and  the  correlation  coefficients  (i>lvpM,pa)  of  the  pooled  data 
are  computed  and  presented  in  Table  8.  furthermore,  the  normalized  pooled  sample  data  in 
Table  7  are  plotted  in  Figure  38  for  vs  C4,  C2  vs  2f4  and  Sf2  vs  S3  to  show  their  correlations. 

It  is  clearly  observed  from  Table  8  and  Figure  38  that  the  correlations  between  C2  and  C3, 
and  between  C2and  C4  are  very  small  whereas  a  strong  correlation  exists  between  C3and  C4  (Bee 
Figure  38(a)).  Consequently,  the  following  assumptions  are  reasonably  made  to  simplify  the 
statistical  analysis  of  C2,  C3  and  C4:  (1)  C2  is  statistically  independent  of  C3  and  C4,  and  (2)  C3 
and  C4  are  completely  correlated  and  hence  they  are  linearly  related: 

Ca  -  «C4  +  0  (43) 

in  which  a  and  0  will  be  determined  later. 

Based  on  these  assumptions,  the  statistical  characterizations  for  C2,  C3  and  C4  become 
tractable,  since  the  establishment  of  a  joint  density  function  for  C2,  C3  and  C4  is  almost 
impossible.  Note  that  the  joint  density  function  is  available  for  C2,  C3and  C4only  when  they  are 
Gaussian  (normal)  random  variables.  It  will  become  apparent  later  that  C2  must  be  positive  and 
hence  it  should  not  be  a  Gaussian  random  variable. 

Talcing  the  expectation  of  Eq.  43  and  forming  the  standard  deviation  of  C3,  one  obtains 

M;,(v,R,T)  -  a/i4(u, R,T)  +  0  (44) 

<r\  (viR.T)  -  <r\  (v,R,T)  (45) 

in  which  «i3(i\R,T)  and  ji4(u,R,T)  are  the  mean  values  of  C3and  C4,  repectively,  and  a3(u,R,T)  and 
(r4(u,R,T)  ore  the  corresponding  standard  deviations. 

The  standard  deviations  in  Eq.  46  can  be  replaced  by  the  mean  value  multiplied  by  the 
coefficient  of  variation  in  order  to  determine  a 

P*  (v,R,T)  V\  -  «*  (t.R.T)  V2  (46) 


M::(v,R,T)V;i 

m4(u,R,T)V4  (47) 


In  Eq.  47,  a  negative  sign  is  chosen  because  the  correlation  between  C3  and  C4  is  negative  (see 
Figure  38(a)). 

Now  0  can  be  obtained  by  substituting  Eq.  47  into  Eq.  44;  with  the  result 

0  -  m:i(v,R,T)  £  1  J  (48) 

It  is  obvious  from  Eqs.  47  and  4S  that  a  *  <*(u,R,T)  and  0  =  0  (v,  R,  T)  depend  on  the  test 
conditions  (v,R,T).  Values  of  a  and  0  under  various  test  conditions  are  computed  from  Eqs.  47 
and  48  as  well  as  Tables  6  and  8,  and  they  are  presented  in  Table  9. 
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TABLE  7.  NORMALIZED  POOLED  DATA  OF  C2  C, 
AND  C4,  M  *=  28 


0.7 8200 E  00 

0.80810E  00 

0.90182E  00 

0.82382E  00 

0.89562E  00 

0.927B6E  00 

0.84263E  00 

0.89682E  00 

0.93891  E  00 

0.86564E  00 

0.91684E  00 

0.94216E  00 

0.87383E  00 

0.92570E  00 

0.94761E  00 

0.89 19BE  00 

0.92882E  00 

0.94B80E  00 

0.93428E  00 

0.95750E  00 

0.96588E  00 

0.94827E  00 

0.9641  IE  00 

0.97226E  00 

0.96842E  00 

0.97936E  00 

0.97422E  00 

0.97084E  00 

0.982 18E  00 

0.97901E  00 

0.97746E  00 

0.98786E  00 

0.98319E  00 

0.98737E  00 

0.99964E  00 

0.98422E  00 

0.99020E  00 

0.10029E  01 

0.98438E  00 

0.99430E  00 

0.10068E  01 

0.98495E  00 

0.99930E  00 

0.10157E  01 

0.99997  E  00 

0.I0006E  01 

0.10212E  01 

0.10028E  01 

0.10090E  01 

0.10226E  01 

0.10030E  01 

0.10122E  01 

0.10246E  01 

0.10046E  01 

0.10162E  01 

0.10245E  01 

0.1014SE  01 

0.I0188F,  01 

0.10267E  01 

0.10176E  01 

0.10645E  01 

0.10387E  01 

0.10214E  01 

0.10853K  01 

0.10477E  01 

0.10287E  01 

0.1 094 IE  01 

0.10699E  01 

0  10416E  01 

0.1109&E  01 

0.10733E  01 

0.10461E  01 

0.11237E  01 

0.10788E  01 

0.10492E  01 

O.I13H6K  01 

0. 10804  E  01 

0.10606E  01 

0.1 2377  E  01 

0.10941E  01 

0.1U16E  01 

0.12637E  01 

0.U330E  01 

0.116S2E  01 

TABLE  8.  COEFFICIENTS  OF  VARIATION  AND 

CORRELATION  COEFFICIENTS  OF  Ca,  C3 
AND  C4  OF  POOLED  DATA 


v, 

Coefficient*  of  Variation 

K, 

V, 

0.1117 

-0,0698 

-0.0545 

Correlation  Coefficient * 

Hi 

Pi, 

0.3384 

-0.9838 

0.3161 
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TABLE  1).  VALUES  OP  «  AND  p  UNDER  VARIOUS 
TEST  CONDITION 


Ti'sl 

Cunditinn 

1 

2 

n 

■I 

,'i 

■If 

■  0, 420H 

-0,4190 

-0,3472 

-0.4433 

-0.4945 

n 

-3.3458 

-3.0972 

-2.844 

-3,3270 

-3,4872 

After  determining  the  mean  values  g2(ii,R,T),  g,,(u,R,T)  and  m4(v,RiT),  and  the  coefficients 
of  variation  (Vj,,V.,,V4)  of  Ca,  C,,and  C4  (Tables  6  and  8),  as  well  as  «  and  f)  (Table  9),  we  are  in 
the  position  to  estimate  the  distribution  functions  for  Cs,  Ca  and  C4.  From  the  physical 
standpoint,  the  crack  growth  rate  da/dn  (or  Y)  should  be  a  monotonically  increasing  function  of 
the  stress  intensity  range  AK  (or  X).  Hence,  it  follows  from  Eq,  1  that  Ca should  be  positive.  We 
shull  select  a  most  appropriate  distribution  function  of  Ca  which  is  defined  in  the  positive 
domain.  Distribution  functions,  such  as  lognormal,  Weibull,  gamma  and  beta,  are  all  defined  in 
the  positive  domuin.  Normalized  pooled  data  presented  in  Table  7  for  C2  are  fitted  by 
(1)  lognormal,  (2)  Weibull,  (3)  gamma,  and  (4)  beta  distributions, 

The  Kolmogorov-Smirnov  (K-S)  test  for  the  goodness-of-fit  indicates  that  the  lognormal 
distribution  fits  the  data  for  Ca  better,  However,  both  the  normal  and  the  lognormal 
distributions  fit  the  data  sets  of  C., and  C4 equally  well.  The  normalized  data  Bet  for  C,4is  shown 
on  the  lognormal  probability  paper  in  Figure  39(a).  The  corresponding  K-S  statistics  Dn  is  also 
given  in  the  figure.  The  normalized  data  sets  for  C,,and  C4  are  plotted  on  the  normal  probability 
paper  in  Figure  U9(b)-B9(c),  along  with  their  corresponding  Dn  values.  The  total  number  of 
pooled  datH  sets  is  28.  With  the  D„  value  given  in  Figure  39,  the  lognormal  distribution  is 
acceptable  for  C,4at  u  significance  level  much  higher  than  20',..  Likewise  the  normal  distribution 
is  acceptable  for  C.,and  C4at  the  same  level  of  significance. 

We  shull  summarize  the  results  obtained  from  the  statistical  analyses  of  the  test  data  for 
U,,and  C4ln  the  following: 

1.  iH  a  lognormal  random  variable.  The  mean  value,  ga(i>,R,T)  as  given  by 
'Fable  (i,  is  a  function  of  the  test  conditions  (u,R,T)  whereas  the  coefficient 
of  variation,  \v  given  in  'Fable  8,  is  independent  of  (i>,R,T). 

2.  U,,  is  statistically  independent  of  C,,  and  C4. 

3.  and  C4  are  normal  random  variables,  The  mean  vulucs,  ^i.t(u,R,T)  and 
M4(m,R,T),  as  given  by  Table  (i,  are  functions  of  t  he  test  conditions  (d,R,T) 
whereas  the  coefficients  of  variation,  V,,  arid  V4,  given  in  Table  8,  are 
independent,  of  (u,R,T). 

4.  C;,  and  C4  are  completely  correlated  and  hence  they  are  functionally  related 
as 

CL,  -  .»(i«, H,T)C4  +  /t(t-,R,T)  (49) 


in  which  «v(i*,R,T)  and  /f(i\R.T)  are  functions  of  test  conditions  (u,R,T)  and  they  are  given  by 
Eijk.  47  and  48,  as  well  uh  'Fable  9. 
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Substituting  Kq.  f>0  into  Kq,  ftf>  and  using  the  theory  of  total  probability,  one  obtains 


Fv(y)  =  fc4  (5)  * 


log  sinh 


■(V)  - 


log  (X  +  +0)  -  m* 


of 


d*  (56) 


in  which  fC4  (£)  is  the  probability  density  function  C4 obtained  from  Eq.  51 

f(:4  (f) 


■r — —  exp'  -  ~  ^  )  ■ 

ir  m4V4  2  V  M4V4  / 

.  • 


(57) 


and  n*j  =  fi*2  ( v.R.T),  <r*2  =  a *s (u,  R,  T)  and  ^4  =  ^(uiR.T)  are  functions  of  (e,R,T)  given  by  Eqs. 
52,  511  and  Table  6.  The  distribution  function  of  the  log  crack  growth  rate  Y  given  by  Eq.  56  can 
be  evaluated  easily  by  a  straightforward  numerical  integration. 

The  distribution  of  the  log  crack  growth  rate,  however,  is  not  the  main  objective  of  our 
research  program.  What  is  important  in  the  present  research  effort  iB  the  life  (or  time)  to  reach  a 
given  crack  Bize.  However,  any  reasonable  statistical  model  should  be  capable  of  describing  the 
physical  behavior  of  test  results.  Hence  nineteen  (19)  sample  functions  of  the  log-crack  growth 
rate  Y  given  by  Eq.  54  have  been  simulated  and  only  the  first  10  of  them  are  displayed  in  Figure 
40.  These  sample  functions  are  simulated  from  Eq.  54  using  the  statistical  distributions  of  C2 
and  C4  under  the  test  condition  No.  1. 

Relation  Batman  Dtatrlbutlona  oi  Crack  Slza  and  Cyelaa  To  Roach  a  Qlvon  Crack  Slza 


The  distribution  of  the  crack  size  a(n,)  at  any  number  of  load  cycles  n,  (or  any  service 
time)  is  the  important  information  needed  in  the  analysis  of  Retirement  for  Cause  for  engine 
components.  Likewise,  the  distribution  of  the  number  of  load  cycles  N(a,)  (or  life)  to  reach  a 
given  crack  size  a,  is  also  very  useful  information.  In  fact,  both  distributions  mentioned  above 
are  related  to  each  other  as  will  be  described  in  the  following, 

^et  ‘«(»|,(x)  be  the  probability  density  function  of  the  crack  size  a(n,)  at  n,  load  cycles,  and 
W,,(n)  be  the  probability  fh-nsity  function  of  the  number  of  cycles  N(a,)  to  reach  a  given  crack 
size  a,.  Both  density  functions  are  schematically  shown  in  Figure  41.  The  probability  that  the 
crack  size  a(n,)  at  n,  load  cycles  will  exceed  a,  is  equal  to  the  probability  that  the  number  of  load 
cycles  N(o,)  to  reach  the  same  crack  size  a,  is  less  than  n,  cycles,  i.e., 


P|a(ni)>ail  =  P(N(a,):<nil 
It  follows  from  Eq.  58  that 

f  f:illl|l(x)dx  -  J"1  fNl„,i(n)dn 


(58) 


(59) 
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Figure  40.  Simulated  Lag  Crack  Growth  Rate  Y  an  Function  of  Lon  Stress  Intensity 
Range  X  Using  Model  2  for  Test  Condition  No.  I 


Equations  58  und  59  indicate  that  the  shaded  areas  shown  in  Figure  41  under  both  density 
functions  are  equal.  Consequently,  if  the  distribution  of  the  number  of  load  cycles  to  reach  a 
crack  size  a,  is  obtained  for  all  values  of  a1(  then  the  distribution  of  the  crack  size  at  any  load 
cycle  n,  (or  service  time)  can  be  derived  from  it,  as  schematically  shown  in  Figure  42.  From 
Figure  42,  the  following  relationships  can  be  used  to  determine  f,(lll)(x)  as  follows: 


J  It  1 

fi,i  It, 

]„'  •’No. 

,  i(n)dn 

(60) 

.c: 

f:.(„li<x)dx 

0  ‘NO. 

,i(n)dn 

(x)dx 

c  <•*., 

,i(n)dn 

82 


n, 

Number  of  Cycles,  n 


ro  3130*3 


Figure  41.  Relation  Between  Probability  for  n,  Load  Cycles  to  Reach  a,  and  Probability 
for  the  Crack  Size  at  n,  To  Be  Larger  Than  at 


Figure  42.  Relation  Between  Distributions  of  Crack  Size  and  Cycles  to  Reach  (liven 
Crack  Sizes 
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Frequently,  the  distribution  function  of  the  crack  size  at  any  service  time  cannot  be  easily 
determined,  but  the  distribution  of  the  time  to  reach  a  crack  size  can  be  estimated  in 
approximation  in  a  relatively  simple  manner.  In  such  a  situation,  the  distribution  of  the  crack 
size  or  the  probability  of  crack  exceedance  will  be  determined  from  that  of  the  time  to  reach  a 
given  crack  size  using  Eq.  60  as  will  be  described  later. 

To  show  the  crack  growth  behavior  based  on  the  present  statistical  model  given  by  Eq.  54, 
referred  to  as  model  2,  nineteen  (19)  simulated  sample  functions  of  the  crack  growth  rate  Y 
presented  above  have  been  integrated  numerically  to  obtain  the  crack  growth  curves  a(n)  versus 
n  using  the  assumed  homogeneous  test  environment  given  in  Table  4  for  the  test  condition 
No.  1.  The  results  are  depicted  in  Figure  43.  It  is  observed  that  the  simulated  sample  functions 
of  the  crack  growth  curves  do  intermingle. 

TEST  CONDITION  NO.  1 


NUMBER  OF  CYCLES/  10**4 

FD  235879 

Figure  43.  Simulated  Crack  Growth  Damage  Accumulation  Using  Model  2  for  Test 
Condition  No.  1 

Statistical  Distribution  of  Time  to  Reach  Given  Crack  Size 

Given  the  statistical  distributions  of  C2  and  C4.  the  distribution  of  the  log  crack  growth  rate 
Y  has  been  derived  analytically  in  Eq.  56.  However,  the  distribution  of  either  the  crack  size, 
a(n),  at  any  number  of  load  cycles,  n,  or  the  number  of  cycles,  N(a),  to  reach  any  crack  size  “a” 
cannot  be  obtained  analytically.  This  is  because  Eq.  54  cannot  be  integrated  in  a  closed  form  to 
yield  the  analytical  relation  among  a(n),  N(a),  C2and  C,.  As  a  result,  the  method  of  perturbation 
will  be  used  to  estimate  the  first  four  ceritralmoments  of  N(a)  as  follows. 
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The  crack  size  a(n)  at  any  number  of  load  cycles,  n,  can  be  obtained  by  numerically 
integrating  Eq.  54  and  expressed  symbolically  as 

a(n)  ---  a(n;Cj,C.|)  (go 


It  is  obvious  that  the  crack  size  a(n;C2,C4)  is  a  statistical  variable,  since  both  C2  and  C4  are 
statistical  variables, 

On  the  other  hand,  Eq.  54  can  be  integrated  numerically  to  arrive  at  the  number  of  load 
cycles  N(a)  to  reach  any  given  crack  size  “a",  which  is  symbolically  denoted  by 


N(a)  -  N(a;C,„C.,)  (62) 

Again,  N(a;C2,C4)  is  a  statistical  variable. 

Let  the  first  four  central  moments  of  C2and  C4be  denoted  by  (M^n22,f2,i72)  and  (M4,a4g,£4,n4)t 
respectively.  Since  C3  and  C4  follow  the  lognormal  and  normal  distributions,  respectively,  their 
first  four  central  moments  can  easily  be  computed  from  the  respective  mean  values  and 
coefficients  of  variation  given  in  Tables  6  and  8.  The  results  are  shown  in  Tables  10  and  11. 
Then,  the  first  four  central  moments  of  N(a),  denoted  by  can  be  obtained 

by  the  method  of  perturbation  as  follows  [Ref.  46], 

MN(a)  =“  N(a;M,,M4)  +  a\ 

"NO  “  (^)  '*2  +  (-lei)  Jff24  +  ($i)(  l63) 

fN<»'  *  ( «&)' 

/ 4  .  /  V  . « /  -<n\2  /  -<n  y ,  , 
m>"  *  W  *  +  h  w  (jcf/  lT* 

in  which  N(a;p2,/i4)  is  the  value  of  N(a)  evaluated  by  integrating  Eq.  54  with  C2*jua,  C4“m4.  and 


<)N  /<N(a;C,,C|) 

<1C,  “  flC- 


£ N  _  «v',N(tt,C.j,C.l) 

«icH =  flCH 


Cj=Mj 

C|*Ml 


iVN  (7N(a;Ca,C.|) 

(lC  |  (7C,| 


C  |“M  l 


fl-  FT  r/'~'N(a;Cj,C.i) 
<ic'i "  <7C- 


C.|=M| 


Since  C4  is  a  normal  random  variable,  the  property  that  has  been  used  in  Eq.  63. 
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TABLE  10.  FIRST  FOUR  CENTRAL  MOMENTS  OF  C2  FOR 
EACH  TEST  CONDITION 


TABLE  11.  FIRST  FOUR  CENTRAL  MOMENTS  OF  C4  FOR 
EACH  TEST  CONDITION 


Tent 

Condition 

S, 

r, 

7, 

1 

-4.47 

o.ofisa 

0 

0.01056 

2 

-4.151 

0.0511 

0 

0.00785 

a 

-4.60 

0.0629 

0 

0.01186 

4 

-4.215 

0.0514 

0 

0.00791 

5 

-a.96 

0.0466 

0 

0.00661 

The  evaluation  of  the  first  fou£  central  moments  for  N(a)  from  Eq.  63  involves  the  partial 
derivatives  ON/dC2),  ON/dC4),  (daN/dC22),  and  (d2N/dCa4)  evaluated  at  C2=M2and  C4=p4,  Eq. 
64.  These  partial  derivatives  can  be  computed  using  the  finite-difference  scheme  as  follows; 


=  lN(a;Mj+A,Mi>  “  N(a;g2-A,fM)l 

~  -  ^  |N(a;M  .,m.1+A)  -  N(a;ju.j,M|-A)|  <66 > 

=  — ,lN(a;ju,+2A.p.|)  -  2N(a;M.>,*<.|)  +  N(a;p2-2A,p.|)l 

dUr;  4  A 

~--7,.[N(a;p.,,p.,+2A)  -  2N(a;a.,M,)  +  N(am2,M.,-2A)l 
(K>“  4  A 


in  which  A  is  chosen  to  be  a  reasonably  small  number  such  that  all  the  partial  derivatives 
converge.  In  the  present  investigation,  the  second  derivative  terms  in  Eq.  63  have  been  neglected 
and  it  iB  found  that  A-0.006  is  sufficient.  Equation  66  involves  various  N(a)  values  which  can  be 
determined  in  a  straightforward  numerical  integration  using  Eq.  64.  For  instance,  the  quantity 
N(a;a2+ A,ji4)  is  estimated  by  setting  C2=fi2+A  and  C4  =  *»4  in  Eq.  64  and  carrying  out  Ihe 
numerical  integration. 

The  first  central  moment  indicates  the  mean  number  of  cycles  to  reach  a  given  crack 
size  “a”,  and  the  second  central  moment  (t2Nui  represents  the  variance  of  N(a).  The  third  central 
moment  fN(<)  is  a  measure  of  skewness  (or  symmetry)  of  the  distribution  whereas  the  fourth 
central  moment  is  related  to  the  peakedness  of  the  distribution  of  N(a),  called  kurtosm.  The 
following  nondimensional  quantities  are  commonly  used  as  measures  of  skewness  and  kurtosis 
of  the  distribution, 
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(66) 


$2  “  (67) 

With  the  assumed  homogeneous  test  environment  given  in  Table  4,  the  mean  values,  j*Nu)1 
the  standard  deviations,  <jN,t),  the  skewness,  \//?7and  the  kurtosis,  (3V  have  been  computed  ana 
presented  in  Tables  12  and  13  for  two  particular  crack  sizes.  The  mean  value  n^t),  the  mean 
value  ±  one  standard  deviation,  gNU)  ±  «rNul,  and  the  coefficient  of  variation,  VN(t)“<rNul/j*N(i,, 
for  test  condition  No.  1  are  plotted  in  Figure  44  and  46,  respectively.  From  these  figures,  aB  well 
as  Tables  12  and  13,  two  interesting  observations  are  made  in  the  following: 

(1)  Although  the  standard  deviation  aN(t)  increases  as  the  number  of  cycles  n 
increases,  the  coefficient  of  variation  VNU)  (a  measure  of  dispersion)  actually 
reduces  aB  n  increases. 

(ii)  The  skewness  \J7TX  of  N(a)  is  quite  small  because  the  perturbation 
approximation  made  in  Eq.  63  takes  into  account  only  the  first  term  in  the 
series  expansion.  Hence  \ffT\  thus  obtained  is  not  very  accurate. 


TABLE  12.  VALUES  OF  DISTRIBUTION  PARAMETER  FOR  WEIBULL,  LOGNOR- 
MAL  AND  GAMMA  DISTRIBUTIONS 


Dintribulinn 

Tent  Crock 

Condition  Sue 

Weibull 

/,<  mnomal 

Gamma 

Perturbationa 

We*  ult* 

«o 

vlnHW 

"inNIul 

7 

\ 

SNi.iI 

"Nllll 

1 

t.c, 

11.788 

69,760 

10.865 

0.289 

11.491 

0.000212 

53,780 

16,018 

■2.(1 

5.102 

75,894 

11.130 

0.219 

20.419 

0.000292 

69,660 

15,447 

2 

0.8 

0.072 

16,990 

9,56 

0.346 

7.902 

0.00062 

16,101 

6,400 

2.0 

4.572 

26,291 

20.056 

0.245 

16.195 

0.00067 

24,000 

6,967 

:i 

0.76 

0.062 

60,972 

10.680 

0.319 

9,291 

0.00020 

46,594 

15,007 

1.7 

6.849 

78,458 

11.100 

0.212 

21.569 

0.00031 

67,620 

14,671 

4 

0.8 

2.612 

18,680 

9.685 

0.300 

10.560 

0.00062 

16,740 

5,171 

2.0 

6.700 

80,420 

10,226 

0.198 

24.938 

0.00086 

28,082 

5,636 

r. 

0.8 

4.704 

9,291 

9.021 

0.235 

17.651 

0.00207 

8,500 

2,026 

2.0 

8,711 

16,498 

9.646 

0.106 

53.279 

0,00341 

15,600 

Theoretically,  the  four-parameter  Johnson  distribution  can  be  employed  for  the  distribu¬ 
tion  of  the  number  of  cycles  N(a)  to  reach  a  given  crack  size  "a.”  The  four  central  moments, 
Mn<«p  ^nui*  an^  7ni*i>  obtained  above  can  be  used  to  determine  the  four  parameters  involved 
in  the  Johnson  distribution.  For  the  practical  analysis  and  design  purpose,  however,  it  is  highly 
desirable  to  simplify  the  approach  as  much  as  possible.  As  a  result,  the  lognormal,  Weibull  and 
gamma  distributions  have  been  chosen  as  possible  candidates  for  the  distribution  of  N(a),  The 
reasons  for  the  choice  of  these  distributions  is  given  as  follows;  (i)  (hey  are  well-known  to 
engineers,  (ii)  they  are  defined  in  the  positive  domain,  and  (iii)  they  involve  only  two 
distribution  parameters. 
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CRACK  SIZE,  A<  N  > 


TABLE  13.  SKEWNESS  AND  KURTOSIS  OF  WEIBULL,  LOGNORMAL  AND 
GAMMA  DISTRIBUTIONS 


Distribution 

Tent  Crack 

Wcibull 

Lognormal 

Gamma 

Perturbations 

Results 

Condition 

Size 

d. 

d> 

d, 

»t 

1 

1.0 

0.041 

2.728 

0.910 

4.510 

0.590 

3.522 

0.126 

3.054 

2.0 

0.280 

2.910 

0.074 

3.820 

0.443 

3.294 

0.124 

3,053 

2 

0.8 

0.14ft 

2.722 

1.112 

5.277 

0.711 

3,76 

0.0548 

3,018 

2.0 

0.190 

2.818 

0.761 

4.046 

0,607 

3.37 

0.0363 

3.010 

a 

0.76 

0.0601 

2.710 

1.019 

4.903 

0.656 

3.645 

0.1369 

3,061 

1.7 

0.209 

2.91)4 

0,8fi6 

3.774 

0.430 

3,278 

0.1560 

3.072 

4 

0.8 

0.002 

2.717 

0.952 

4.653 

0.015 

3.567 

0.0738 

3.027 

2.0 

0.360 

,1.00!) 

0.608 

3.666 

0.400 

3.241 

0.0517 

3.016 

n 

(1.8 

0.224 

2,850 

0.727 

3.955 

0,476 

3.340 

0.110 

3.046 

2.0 

0.R74 

8.417 

0.418 

3.306 

0,274 

3.112 

0.086 

3.032 
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Figure  44,  Average  Number  of  Load  Cycles  uNl0,  to  Reach  a  Crack  Size  and  Average 
One  Standard  Deviation  nN(aj  ±  aNlul  for  Test  Condition  No.  1 
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TE8T  CONDITION  NO.  1 
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Since  each  distribution  function  lifted  ahove  involves  only  two  parameters,  the  first  two 
central  moments  of  N(a),  #iNll)1and  <TvlN(lt),  obtained  from  the  method  of  perturbation  (eq.  fill)  are 
used  to  estimate  the  two  distribution  parameters  using  the  method  of  moment.  The  distribution 
parameters  are  shown  in  Table  12  for  each  test  condition  (i  .ll.T)  at  two  crack  sizes. 

After  determining  the  distribution  parameters  for  each  distribution  function  using  the  first 
two  central  moments  from  the  perturbation  results,  the  skewness  y/js,  and  the  kurtosis  f^can  be 
computed  from  the  respective  distribution  parameters.  They  are  presented  in  Table  13  along 
with  the  perturbation  results.  Then,  the  goodness-of-fil  for  each  distribution  may  be  judged  by 
comparing  the  v/tf,  and  values  of  each  distribution  with  those  of  the  perturbation  results.  It  is 
observed  from  Table  13  that  the  Weibull  distribution  is  the  best  of  the  three  distributions 
considered  above.  Thus,  the  Weibull  distribution  is  selected  for  the  distribution  of  N(a),  the 
number  of  cycles  to  reach  a  given  crack  Bize  "a." 

Condition  With  Tttl  R*»ult» 

A  qualitative  correlation  study  Is  made  of  the  second  statistical  model  with  respect  to  the 
extrapolated  test  results  in  this  section.  Although  the  statistical  distribution  of  the  log  crack 
growth  rate  Y  hus  been  derived  in  Eqs,  56  and  57,  it  would  be  more  relevant  to  examine  the 
behavior  of  the  sample  functions  of  Y.  In  this  connection,  19  sample  functions  of  Y  versus  X 
have  been  simulated  using  Eq.  54  with  the  statistics  of  C,  (1= 2,3,4)  given  by  Tables  6  and  B,  The 
results  are  presented  In  Figures  46-50  for  all  test  conditions.  With  the  assumed  homogeneous 
test,  environments  given  in  Table  4,  the  craok  growth  rates  displayed  in  Figures  46-50  have  been 
integrated  to  yield  the  simulated  oruck  growth  damage  accumulation.  The  corresponding  crack 
size  u(n)  versus  the  number  of  load  cycles  n  are  shown  in  Figures  51-55. 

The  average  number  of  load  cycles,  gN(lll,  to  reach  any  given  crack  Hize  and  the  standard 
deviation,  <tNU),  have  been  computed  using  the  method  of  perturbation.  The  results  for  mN(11)  and 
^N(»i  ^nui  «re  shown  in  Figures  56-60.  A  comparison  between  the  corresponding  results  in 
Figures  56-60  and  Figures  lH(h),  23(b)-26(b)  indicates  that  the  correlation  between  the  second 
statistical  model  and  the  extrapolated  test  results  is  reasonable, 

Ah  described  previously,  the  Weibull  distribution  is  suitable  for  the  distribution  of  the 
number  of  load  cycles  to  reach  any  crack  size.  With  the  Weibull  parameters  nr,, and  d0  estimated 
in  Table  12,  the  dist  ributions  of  N(a)  for  two  cruck  sizes  is  depicted  as  dashed  curves  in  Figures 
27-31.  Likewise,  the  crack  exceedance  curves  based  on  the  present  model  are  shown  in  Figures 
32-36  uh  dashed  curves.  It  is  observed  from  these  figures,  that  the  correlation  between  the 
second  statistical  model  (dashed  curves)  and  the  extrapolated  test  results  (circles)  is  also 
reasonable. 
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X-  LOG< DELTA  K> 
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Figure  46.  Simulated  Log  Crack  Growth  Rate  Y  as  Function  of  Log  Stress  Intensity 
Range  X  Using  Model  2  for  Test  Condition  No.  1 
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Figure  47.  Simulated  Log  Crack  Growth  Rate  Y  as  Function  of  Log  Stress  Intensity 
Range  X  Using  Model  2  for  Test  Condition  No.  2 
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Figure  50.  Simulated  Log  Crack  Growth  Rate  Y  as  Function  of  Log  Stress  Intensity 
Range  X  Using  Model  2  for  Test  Condition  No.  5 
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Figure  51.  Simulated  Crack  Growth  Damage  Accumulation  Using  Model  2  for  Test 
Condition  No.  I 
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figure  52.  Simulated  Crack  Growth  Damage  Accumulation  Using  Model  2  for  Test 
Condition  No.  2 
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Figure  54.  Simulated  Crack  Growth  Damage  Accumulation  Using  Model  2  for  Test 
Condition  No.  4 
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Figure  55.  Simulated  Crack  Growth  Damage  Accumulation  Using  Model  2  for  Test 
Condition  No.  5 
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TEST  CONDITION  NO.  1 
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Figure  58.  Average  Number  of  Load  Cycles  fiN(a)  to  Reach  a  Given  Crack  Size  and 

Average  ±  One  Standard  Deviation  gNM  ±  aN(t)  for  Test  Condition  No.  3 
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Figure  59.  Average  Number  of  Load  Cycles  uN,al  to  Reach  a  Given  Crack  Size  and 

Average  ±  One  Standard  Deviation  aN(tl)  ±  „N(a)  [or  Test  Condition  No.  4 
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Figure  60.  Average  Number  of  Load  Cycles  iiN(a)  to  Reach  a  Given  Crack  Size  and 

Average  ±  One  Standard  Deviation  uNlaJ  ±  «NW  for  Test  Condition  No.  5 
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Since  each  distribution  function  listed  above  involves  only  two  parameters,  the  first  two 
central  moments  of  N(n),  /jN(l0and  <t'2n<<1,,  obtained  from  the  method  of  perturbation  (eq.  63)  are 
used  to  estimate  the  two  distribution  parameters  using  the  method  of  moment.  The  distribution 
parameters  are  shown  in  Table  12  for  each  test,  condition  (i>,R,T)  at  two  crack  sizes. 

After  determining  the  distribution  parameters  for  each  distribution  function  using  the  first 
two  central  moments  from  the  perturbation  results,  the  skewness  \JJS{  and  the  kurtosis  /32  can  be 
computed  from  the  respective  distribution  parameters.  They  are  presented  in  Table  13  along 
with  the  perturbation  results.  Then,  the  goodness-of-fit  for  each  distribution  may  be  judged  by 
comparing  the  \J(3 ,  and  fi2  values  of  each  distribution  with  those  of  the  perturbation  results.  It  is 
observed  from  Table  13  that  the  Weibull  distribution  is  the  best  of  the  three  distributions 
considered  above.  Thus,  the  Weibull  distribution  is  selected  for  the  distribution  of  N(a),  the 
number  of  cycles  to  reach  a  given  crack  size  “a.” 

Correlation  With  Test  Results 

A  qualitative  correlation  study  is  made  of  the  second  statistical  model  with  respect  to  the 
extrapolated  test  results  in  this  section.  Although  the  statistical  distribution  of  the  log  crack 
growth  rate  Y  has  been  derived  in  Eqs.  56  and  57,  it  would  be  more  relevant  to  examine  the 
behavior  of  the  sample  functions  of  Y.  In  this  connection,  19  sample  functions  of  Y  versus  X 
have  been  simulated  using  Eq.  54  with  the  statistics  of  C|  (i=2,3,4)  given  by  Tables  6  and  8.  The 
results  are  presented  in  Figures  46-50  for  all  test  conditions.  With  the  assumed  homogeneous 
test  environments  given  in  Table  4,  the  crack  growth  rates  displayed  in  Figures  46-50  have  been 
integrated  to  yield  the  simulated  crack  growth  damage  accumulation.  The  corresponding  crack 
size  a(n)  versus  the  number  of  load  cycles  n  are  shown  in  Figures  51-55. 

The  average  number  of  load  cycles,  /xN(a),  to  reach  any  given  crack  size  and  the  standard 
deviation,  <xN(a),  have  been  computed  using  the  method  of  perturbation.  The  results  for  /iN(a)  and 
MN(a)  ±  <rN(a)  are  shown  in  Figures  56-60.  A  comparison  between  the  corresponding  results  in 
Figures  56-60  and  Figures  18(b),  23(b)-26(b)  indicates  that  the  correlation  between  the  second 
statistical  model  and  the  extrapolated  test  results  is  reasonable. 

As  described  previously,  the  Weibull  distribution  is  suitable  for  the  distribution  of  the 
number  of  load  cycles  to  reach  any  crack  size.  With  the  Weibull  parameters  n^and  /?0  estimated 
in  Table  12,  the  distributions  of  N(a)  for  two  crack  sizes  is  depicted  as  dashed  curves  in  Figures 
27-31.  Likewise,  the  crack  exceedance  curves  based  on  the  present  model  are  shown  in  Figures 
32-36  as  dashed  curves.  It  is  observed  from  these  figures,  that  the  correlation  between  the 
second  statistical  model  (dashed  curves)  and  the  extrapolated  test  results  (circles)  is  also 
reasonable. 
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